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ABSTRACT 

We have developed antenna-coupled transition-edge sensor (TES) bolometers for a wide range of cosmic mi¬ 
crowave background (CMB) polarimetry experiments, including Bicep2, Keck Array, and the balloon borne 
Spider. These detectors have reached maturity and this paper reports on their design principles, overall per¬ 
formance, and key challenges associated with design and production. Our detector arrays repeatedly produce 
spectral bands with 20%-30% bandwidth at 95, 150, or 220 GHz. The integrated antenna arrays synthesize 
symmetric co-aligned beams with controlled side-lobe levels. Cross-polarized response on boresight is typi¬ 
cally ~ 0.5%, consistent with cross-talk in our multiplexed readout system. End-to-end optical efficiencies in 
our cameras are routinely 35% or higher, with per detector sensitivities of NET^300 /iKcMBv^. Thanks to the 
scalability of this design, we have deployed 2560 detectors as 1280 matched pairs in Keck Array with a com¬ 
bined instantaneous sensitivity of ^ 9 /iKcMB\/s, as measured directly from CMB maps in the 2013 season. 

Similar arrays have recently fiown in the Spider instrument, and development of this technology is ongoing. 

Subject headings: cosmic background radiation — cosmology: observations — instrumentation: polarime- 
ters —detectors: antenna coupled TES bolometers 
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Cosmic microwave background (CMB) polarimetry is a key 
observable to further our understanding of cosmology in both 
the later and early universe. Degree-scale B-mode polariza¬ 
tion can be used to constrain the tensor-scalar ratio r and 
plac e limits on the energy scale and potential form of infia - 

tion ( [Zaldarriaga & Seljak|1997[[Kamionkowski et al.|1997| ). 

Arcniinute B-mode measurements allow precise reconstruc- 
tion of the gravitational leasing potential at later times, which 
in turn can constr ain the neutrino masses a nd the dark energy 
equation of state ( [Kaplinghat et al.||2003| ). Precise measure¬ 
ments of E-mode polarization provide further cosmological 
information on the plasma physics at recombination. Very 
precise future lensing polarization maps could ultimately be 
used for deeper searches of infiationary polarization. Einally 
a future space mission has the potential to measure large-scale 
polarization to astrophysical limits, for precise tests of infia- 
tion. However, these ambitious scientific goals require high 
sensitivity cameras with exquisite control of systematic er¬ 
rors. 

To meet this need, we have developed large arrays of 
dual-polarized antenna-coupled transition edge sensor (TES) 
bolometers. The key feature in our design is optical coupling 
through a planar antenna, allowing the entire design to be fab¬ 
ricated with scalable photolithographic techniques. This has 
allowed us to rapidly deploy arrays for Bicep2, Keck Array, 
and the balloon borne Spider, amounting to over 6000 de¬ 
tectors fielded as of this writing. We deployed Bicep2 for 
observing in 2010, three Keck Array 150 GHz cameras for 
observing in 2011, and five 150 GHz cameras for Keck Ar¬ 
ray in 2012. In 2014, we replaced two Keck Array 150 GHz 
cameras with 95 GHz cameras, and in 2015 we replaced two 
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more with 220 GHz cameras (fielding an additional 1000 de¬ 
tectors). Spider recently conducted a long-duration flight 
from Mc Murdo Station in ear ly, deploying an additional 2400 
detectors ( [Fraisse et al.|2013 ). 

This paper describes the design principles of these detec¬ 
tors, as well as their on-sky performance and some challenges 
associated with their production. We draw from our exten¬ 
sive experience with the aforementioned experiments, and de¬ 
scribe the as-deployed performance of our devices. 

The format of this paper is as follows. In ^ and ^ we 
describe the cameras and detector element designs, ^bo- 
rating beyond the material in t he Bicep2 Instrument Paper 
( [Bicep2 Collaboration II|2014] ). ^describes the fabrication 
techniques used. To optimize our designs and recipes, we had 
to characterize the optical properties of our thin films for mil¬ 
limeter waves, which we describe in ^ ^describes the an¬ 
tenna array beam synthesis as well as related processing and 
design challenges. Our high on-sky detector yield is only pos¬ 
sible through control of detector properties across each tile, 
which we discuss in ^ ^describes the attained camera sen¬ 
sitivity and we show rough agreement between those mea¬ 
surements and a noise model. Finally, ^ offers some con¬ 
cluding remarks and describes future endeavors. 

2. OVERVIEW 

Cameras — Because the detector arrays are coupled to ex¬ 
periments with refracting optics, we briefly summarize their 
architecture. The Bicep2, Keck Array, and Spider cam¬ 
eras are similar f/2.2 26 cm aperture telecentric refracting 
telescopes, which re-image the sky onto the focal plane with 
pairs of high-density polyethylene (HDPE) lenses. Bicep2 
and Keck Array contain teflon, nylon, and reflective metal- 
mesh edge-filters that filter incident radiation to maintain low 
thermal loading on the focal plane such that internal closed- 
cycle ^He/^ He/'^He sorption fridges can cool the focal plane 
to 280 mK dPuband & Collaudin|1999| ). Due to loading con¬ 
straints at float, the Spider cameras cannot employ lossy 
teflon filters and instead utilize metal-mesh filters. They still 
use a nylon filter on the cold side of 4K. Spider uses a ^ 2 K 
pumped helium bath and single stage closed cycle ^He fridges 
to chill its focal planes to ^ 300 mK. Spider uses ultra-high- 
molecular-weight polyethylene (UHMWPE) windows in its 
low pressure flight environment, whereas Bicep2 and Keck 
Array use closed cell nitrogen filled Zotefoam (respectively 
Propozote PPA30 and Plastazote HD30). We refer the reader 
to the Instrument Paper for m ore BlCEP2-specific details •HE 
|cep2 Collaboration II|2014| ). 

Focal plane — Each focal plane (Figure contains four de¬ 
tector tiles mounted to a common copper frame. Each tile is 
a square of silicon cut from a standard 100 mm wafer, typ¬ 
ically 550 yum thick. The 150 GHz Bicep2, Spider, and 
Keck Array focal planes contain four tiles, each consisting of 
64 dual orthogonally polarized detector pairs, providing 512 
detectors per camera. The detectors’ beams terminate on their 
camera’s cold stop at -15 dB of the main lobe with 1.8/A 
spacing. The 95 GHz Keck Array cameras have tiles con¬ 
taining 36 orthogonally polarized detector pairs, and thus 288 
detectors per camera. These detectors’ beams terminate on 
their camera’s stop at -12 dB of the main lobe with 1.5/A 
spacing. Current 220 GHz focal planes contain four tiles of 
64 detector-pairs each in order to match a readout system de¬ 
signed for the 150 GHz cameras; this under samples the focal 
plane (2.6/A spacing), which will be addressed with a future 


higher density focal plane design. The 220 GHz detectors’ 
beams are similar in size to the 150 GHz beams, terminating 
on the stop at -15 dB. 

We mount the tiles to the aforementioned gold-plated cop¬ 
per frame and thermally sink the tiles with gold wire bonds 
linking the frame to gold bond-pads that directly contact the 
silicon substrate of the tile. The antennas are patterned on the 
non-illuminated side of the wafer, so th at lig ht arrives at the 
antenna through its silicon substrate (see |3.1| for elaboration). 
A quartz anti-reflection tile is mounted on the illuminated side 
of the wafer, while the antenna side faces a superconducting 
niobium reflective back-short placed A/4 away. Testing of 
early Bicep2 prototype focal planes suggested that near the 
tile edge, antennas for one of the two polarizations may cou¬ 
ple to the frame, yielding elliptical beams in the far field. We 
have mitigated this coupling in all deployed focal planes by 
keeping the edge detector pairs a few wavelengths away from 
the frame and by corrugating the frames with specifically cho¬ 
sen depth and impedance grooves, optim ized with sim ulations 
performed in the CST Microwave Suite ( [Corp^OO^al ). 

Polarized Detector Elements — Our detector design is entirely 
planar and does not require horns or other contacting optics. 
In each detector element, optical power couples to two co¬ 
located, orthogonally polarized planar antenna arrays, each 
composed of slot sub-radiators patterned in a superconduct¬ 
ing niobium (Nb) ground plane. All slots of a given orienta¬ 
tion are coherently combined through a micro strip summing 
tree to synthesize a single equivalent antenna for that polariza¬ 
tion orientation. Power from each antenna is passed through 
an on-chip band-defining filter before being dissipated on a 
suspended bolometer island. A superconducting TES on that 
island detects variations in the power received by the anten¬ 
nas. 

3. DETECTOR DESIGN 
3.1. Antenna Design 

The antenna slots in each detector must be spaced to 
Nyquist sample the focal plane surface to avoid grating lobes 
that would rapidl y change the impedance with frequency 
( |Kuo et aL]|2008| ). The antenna pattern of each axis of an 
array is calculated from the N elements per linear dimension 
spaced at distance 5' as follows: 

(V-l)/2 

A(6»)= Y, 

m=-{N-\)/2 

^m{N^TS^Je'r^mO /\o) 

/\o) 

where A^ is the free-space wavelength, the relative permit¬ 
tivity of the surrounding medium, and the sum is across sub¬ 
antennas indexed by m . In addition to the strong peak in the 
normal direction (6> = 0), there are grating lobe peaks when 
^Je'rS^m{0)/\o is a positive integer. To avoid these lobes, the 
slot spacing must be 



where Xo,mm is the minimum wavelength of operation and 
the term in parentheses accounts for the finite width of the 
grating-lobe peaks. For the 150 GHz detectors fabricated 
on silicon (e^ = 11.8) with an upper band edge of 180 GHz 
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Eig. 1.— Eocal Plane pictures and cross-section. Upper. Top (left) of the Bicep2 focal plane and bottom with backshort removed (right). The arrays of 64 
detector pairs per tile are visible at top right. Lower. Major component layers of the focal plane design, with an expanded view of the tile layers at right. Gold 
wire-bonds thermally sink the tiles to the frame and the frame bears corrugations to suppress coupling to the detectors. Only the detectors at the tile perimeter are 
at risk of frame coupling. Corrugations are visible in the inset photo in the upper left panel. 


(K,min = 1.7 mm), the spacing must satisfy 5- < 460 /im. To 
achieve this high density of radiators without intersection, we 
couple power to offset (echelon) pairs of radiators with a Bra- 
vais lattice defined in Figure These dual-slot sub-arrays 
tile the detector elements in 8x8, lOx 10, or 12x 12 ver¬ 
sions, where the overall detector size is chosen to match the 
f/2.2 camera optics. By using such a large number of sub¬ 
antennas, we avoid excessive excitation of substrate modes 
that might degrade the detector’s efficiency. The offset slot- 
pair geometry allows the two orthogonally polarized antenna 
arrays within each detector pair to be co-located (visible in 
upper two panels of Figure®. 

It is energetically favorS)le for the antennas to receive 
power through the silicon substrate to the vacuum side at a 
ratio of 377 O : |Z^”|, where (in equation]^ differs from 
silicon’s TEM impedance of 110 Q in phase because standing 
waves in the substrate modify the effective impedance seen 
by the slots. We exploit this power difference by orienting the 
tiles with the silicon substrate side towards the sky and then 
terminating the back response on the vacuum side with a A/4 
back-short. We mount A/4 quartz anti-refiection (AR) tiles to 
minimize refiection at the air-substrate interface. 

3.2. Antenna Impedance 

Optical power couples from the slots to a microstrip feed 
network. To minimize return loss at this interface, we need 
to accurately compute the input impedance of that feed net¬ 
work. We have written custom software that provides in¬ 
creased speed and versatility for this task over commercial 


options. 

We compute the input impedance of our antenna from the 
required jump-discontinuity in magnetic fields across the slot 
at the microstrip feed-points (/feed = where we impose a 
driving current Jy. Modeling the currents in the ground plane 
with fictitious magnetic currents M = x E running longitu¬ 
dinally down the slot, the field jump-discontinuity requires: 

Jy = f dx'dy'Mix',y') 

[Ghu,M x(.x-x',y-y')] 

The antenna impedance can be inferred once these currents 
M are solved for. The Green functions describe the 

fields on either side of the slot (^ = [1,2]]) resulting from an 
infinitesimal magnetic current in the x direction (parallel to 
the slot.) In the spectral domain, these have the form 


GH,MSkx^ky,Z=0)-- 


cos((l)k) ^ sin((t)k) 

4m(^p) 4e(^p) 


(4) 


where (pk = 2ivci2in{ky/kx) and kp = ^k^+k^ (Das & Pozar 

|1987| ). The dielectric films support transverse electric 
(TE) and transverse magnetic (TM) modes with different 
impedances. The impedance of each mode “seen” by 

the slots can be computed by solving for the standing waves 
of electric and magnetic fields in a manner analogous to treat¬ 
ing the silicon substrate and quartz AR coating as transmis¬ 
sion lines and then transf orming the impedan ce of free space 
(377 Q) through them. ( |Das & Pozar||1987] ) The backshort 
(ground) transforms through a quarter-wave space, presenting 
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Fig. 2.— Optical and SEM photographs of major features of the antenna 
array. Upper, one quarter of a detector element. The antenna array, one filter, 
and one TES bolometer are visible, as well as DC readout lines for the detec¬ 
tor. Middle: SEM micrograph of the slot array (dark rectangles) and oblique 
Bravais lattice (arrows). The thin white lines comprise the microstrip feed. 
For 150 GHz detectors, a~ 600 pm and h = a/2 ^ 300 pm\ the slot dimen¬ 
sions and spacing in the 95 GHz are 63% larger and those in the 220 GHz 
detector elements are 47% smaller. Lower: SEM micrograph of microstrip 
crossover and shunt capacitor at a sub-antenna slot. 


a parallel impedance that is open at the band center. The sili¬ 
con substrates are a standard 100 mm diameter 550 fim thick 


wafer, and we optimized the 95 GHz and 150 GHz impedance 
matching with this thickness in mind. The 220 GHz versions 
are fabricated on thinner 370 /im thick silicon substrates to 
scale this design. 

We expand the unknown fictitious currents M(x,y) = 
y^^iVjFiix^y) in a one-d imensional piecewise sinusoidal basis 
iKominaini et al.|1985] ): 


Fi(x,y) = 


7rA/(w/2)2-3;2 1 0, 


sin(kp(Ax-\x-Xi \)) 
sin(kp Ax) ’ 


Otherwise 


< Ax 


(5) 

localized in physical space at position Xi over the slot, with 
sections of length Av and width w corresponding to the slot 
width. This basis models the expected current distribution on 
the shorter slot dimensions, leaving the unknown distribution 
only on the longer dimension. Moreover the Fourier trans¬ 
form of the basis function has a compact analytic form. These 
features greatly simplify the computations. Solving Equa¬ 
tion [3] with a Galerkin moment method reduce s it to a matrix 
equation K 7 ][V;] = [Ij] ( Kominami et al.|1985 ), where 




dk,dkyF;(k,,ky)G(k,,ky,Z = 0)Fj(k,,ky) ( 6 ) 


The inverse of the admittance matrix in Equation is an 
impedance between basis functions. The self and mutual 
impedance between basis functions located at feed-points is 
the input impedance that our feed-lines must match. 

We avoid complicated numerical integration and further re¬ 
duce the number of unknown V^s by assuming the arrays are 
infinite in extent and applying periodic boundary conditions 
where the antenna properties repeat for spectral translations 
along the lattice vectors: 


where a and b are the distances between the slots in the hori¬ 
zontal and vertical axes (see Figure [^. This approximation is 
valid for slots in the interior of the array and reduces the inte¬ 
grals to a discrete sum that we terminate at a wavevector high 
enough to allow convergence. Figure shows the computed 
input impedance against frequency for the 150 GHz antennas. 
The oscillatory behavior at 200 GHz is expected and real {i.e. 
not a numerical artifact) and coincides with the onset of grat¬ 
ing lobes. We note that we have performed more numerically 
intensive simulations of finite slot arrays where continuous in¬ 
tegration retains edge effects to some degree and these agree 
well with the infinite array approximation. 

We couple power to microstrip lines at the points near the 
slot ends where the lines cross and shunt to the ground-plane 
on the opposite side (see lower panel of Figure [^. This 
coupling effectively transmutes the electric fields across the 
slot into fields between the microstrip conductor and ground 
plane, provided that impedances are well matched. To avoid 
grating lobes, we keep the slots as short as possible and excite 
the slots at their first resonance. However, the radiation re¬ 
sistance of center-fed slots at this resonance is nearly 300 
caused by a current node at the slot center. By symmetri¬ 
cally feeding the slots with a pair of off-center microstrip lines 
close to the current anti-nodes, we greatly reduce the radi¬ 
ation resistance to ^40 Ft. The microstrip feed lines must 
match this real antenna resistance, which, for our dielectric 
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and metal film thicknesses, corresponds to the thinnest lines 
that we can reproducibly fabricate. The feed points also have 
a 5 inductive reactance, which is relatively stable over the 
band. We tune this away with series capacitors that shunt cur¬ 
rent to ground (also visible in Figure |^. Future designs may 
use more complicated terminations {e.g ., open stubs) in lieu of 
capacitors to expand the bandwidth, allowing for co-located 
dual-color detectors. 



Fig. 3.— Simulated feed-point antenna impedance vs. frequency. The 
dashed vertical line is the band center. 


3.3. Micro strip feed 

The waves from the slots coherently sum in the microstrip 
feed-network, which accomplishes beam synthesis in lieu of 
a horn. Figure summarizes this topology. This network 
combines waves across rows, and then sums the waves from 
each row in a column tree at the side of the detector element, 
presenting a single microstrip line in each polarization. The 
two orientations of slots couple power to two independent 
feeds that tightly interlock across the detector pair. The feeds 
are meant to be corporate, combining all waves with uniform 
phase, but the trees are not necessarily binary. In fact, the 
220 GHz and 150 GHz detectors have a 12 x 12 array for¬ 
mat while the 95 GHz ones are 10 x 10, neither of which is a 
power of 2. 



Fig. 4.— Abbreviated feed-network schematic for one detector pair, show¬ 
ing the branch structure for the two polarization summing trees. 

Waves sum in microstrip tee-junctions, with impedances 
chosen to match across each junction when looking from the 


port closest to the bolometer. We pick the ratio of impedances 
on the ports closest to the slots to determine the illumina¬ 
tion pattern. To accomplish this synthesis with low return 
loss, we must construct microstrips with correct impedances, 
which requires accurate knowledge of dielectric constants and 
the penetration depth into our superconducting niobium films. 
(See Thus far, all detectors deployed for Bicep2 Keck 
Array and Spider produce a uniform top-hat illumination 
where the power splits are in proportion to the number of 
slots on each side of the junctions. The microstrip impedances 
can be chosen to synthesize an arbitrary pattern, however, and 
Gaussian tapered feeds have recently been deployed in the 
Biceps camera. 

Coupling between microstrip lines can generate phase er¬ 
rors across the antennas, which we elaborate on in ^ The 
present designs have judiciously spaced the adjacent lines to 
minimize this coupling. To negate any residual phase errors, 
the feed-networks also contain sections of transmission line 
before each slot whose length we can vary during fabrication. 

3.4. Band-defining filters 

Each microstrip feed contains an integrated band-defining 
filter between the antenna feed and bolometer. We use a three- 
pole design as a compromise between bandwidth and loss: 
additional poles would let us increase the bandwidth and still 
avoid atmospheric features, but the additional passes of the 
waves through each resonator increase filter loss. Figure 
shows a micrograph of the filter and a schematic of its equiv¬ 
alent circuit. We realize the resonators with lumped compo¬ 
nents, which do not suffer from the high frequency resonant 
leaks present in A/2 and A/4 transmission line resonators. 



Fig. 5.— Microscope photograph of hlter and equivalent circuit for 
150 GHz. 


Each pole is a series LC resonator, in which the induc¬ 
tors are short stretches of high-impedance coplanar waveg- 
uide (CPW). The CPW impedance Z = .JlJc of roughly 50 0 
exceeds the surrounding lines and acts as a series inductor 
by allowing strong magnetic fields in the CPW ground-plane 
gaps. This inductance is almost entirely magnetic; kinetic in¬ 
ductance from the niobium superconducting microstrip and 
ground layers only makes a minor correction. The series ca¬ 
pacitors are parallel-plate metal-insulator-metal between up¬ 
per and lower niobium films using the microstrip Si 02 as 
the dielectric. In an ideal design the central resonator would 
be a parallel LC resonator shunting to ground. For ease of 
fabrication we instead build a series resonator and invert its 
impedance on resonance using shunt capacitors to gr ound and 
reduced series capacitors (jGalbraith & Rebeiz|2008|). 

We numerically optimize the filter in si mulations with 
the commercially available Sonnet software ( [Corp|[^09b| ). 
The optimizations are constrained to maintain an entrance 
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impedance of 10 which results in fabricable capacitors and 
inductors. Our design results in inductors longer than A/8 that 
have some slight parasitic shunt capacitance, but this does not 
degrade the performance appreciably. The filter in Figure 
is tuned for a 150 GHz band-center. Filters for 95 GHz ana 
220 GHz versions have components with the same reactance 
at the resonant frequencies. We show spectra and summarize 
spectral response features in ^ 

3.5. Bolometer Design 

We terminate and thermalize millimeter-wave power on 
a released bolometer island in a meandered lossy gold mi¬ 
crostrip. 150 GHz waves are constrained to propagate through 
a skin depth of 180 nm, and over a length of several wave¬ 
lengths the line absorbs 99% of the power between two 
passes. We scale the 220 GHz and 95 GHz designs appro¬ 
priately to maintain -20 dB return loss. 


Silicon-Nitride Legs 



Fig. 6.— Electron micrograph of a released TES bolometer, illustrating its 
major components. The gold-meandered microstrip termination is at the right 
of the photograph and the TESs at left. The thicker gold him in the center of 
the island ensures thermal stability. 


This termination is in close thermal contact with two tran¬ 
sition edge sensors: an aluminum TES with a transition tem¬ 
perature Tc ^ 1.2 K for lab tests and a ^ 60 mfl titanium TES 
with a Tc ^ 0.5 K for on-sky observations. The saturation 
power (optical plus electrical power) needed to bring the TES 
temperature to Tc is 


Psat = GrTr 


n-\-l 


( 8 ) 


where the conductance Gc is evaluated at the transition tem¬ 
perature Tc and the surrounding heat bath is Tg ^ 280 mK. 
The exponent n reflects the thermal carriers in the legs, where 
n=\ would correspond to electrons and ^ = 3 would corre¬ 
spond to phonons in 3-D materials. Our typical devices are 
described by ^ = 2.5. The aluminum Tc affords the bolometer 
a high saturation power for use under a 300 K background, but 
has higher noise that would be unacceptable for astrophysical 
observations. 

Each bolometer island is suspended from the tile by six iso¬ 
lation legs: one carrying the microstrip from the antenna, one 
carrying the TES DC lines, and four thinner legs for mechani¬ 
cal stability. Some newer devices use four legs, while Spider 
employs long meandered legs as described below. We tune 
the bolometer legs to achieve the design thermal conductance. 
We aim to keep the phonon noise subdominant to the photon 
noise, but also to keep the G high enough that the detectors do 


not saturate under typical on-sky loading conditions, which 
for Bicep2 is 4-6 pW. We typically build in a safety factor 
of 2 beyond the expected maximum loading; for Bicep2 this 
resulted in Gc = 80-150 pW/K, as described in ^ Two de¬ 
tector pairs in the corners have the microstrip opened between 
the antenna and bolometers such that the TESs do not receive 
opt ical stimulation through the microstrip feed. As described 
in |6.3| we use these four bolometers for diagnostic measure¬ 
ments and to monitor direct stimulation of our detectors by 
optical power that bypasses the microstrip circuitry and may 
thus not have the required spectral and polarization properties. 

Thick evaporated gold (visible in Eigure is added to 
the bolometer islands to boost the heat-capacity by Cau — 
0.5 -0.3 pJ/K, bringing the time constant r = C/G to around 
a millisecond. We infer this from responses to a 1 Hz square- 
wave modulated broad-spectrum noise signal. The additional 
gold is important because of the strong electrothermal feed¬ 
back of these devices; without it, our detectors would enter 
electrothermal oscillations throughout much of the TES tran¬ 
sition. 

We voltage bias the TESs into the transitions typically at 
half the normal resistance {i.e. /?bias — ^n/2), using 3 mD 
parallel shunt resistors. We detect changes in the TES current 
using a SQUID-based time-domain multiplexing architecture, 
where all detectors in a set of 32 are sequentia lly r ead through 
a comm on set of lines ( |de Korte et al.| ( |200^ and |Stiehl et al.| 
( |2011b ). Eor Bicep2, our revisit rate of a given detector is 

25 kHz, ensuring that we Nyquist sample the time-stream 
up to the 6 kHz roll-off imposed by a 1.35 /iH series inductor. 

Adapting this technology to the low photon loadings of Spi¬ 
der’s long-duration balloon payload poses challenges in de¬ 
tector and instrument design. In order to take full advantage 
of the low photon noise levels available at 36 km altitudes, 
care must be taken to ensure that the bolometers achieve low 
noise levels and the surrounding instrument contributes mini¬ 
mal additional photon loading. 

Spider employs several modifications to the Bicep2 de¬ 
tector design to reduce detector noise. A meandered leg de¬ 
sign reduces conductance to G ^ 15 pW/K, reducing phonon 
noise while remaining within a similar footprint. This chosen 
G is not limited by the anticipated loading of ^ 4.5 Kcmb but 
rather to have margin on the 300K background in the labora¬ 
tory while biased on the Al TES. Lower G leads to naturally 
slower detectors, reducing the need to add heat capacity to 
the island with thick gold (Au) to maintain stability. The long 
cabling in Spider prevents us from multiplexing faster than 
^ 20 kHz, so we control aliasing with larger inductors. 

4. FABRICATION 

We fabricate detectors in monolithic batches of 64 matched 
pairs at the Microdevices Laboratory (MDL) at the Jet Propul¬ 
sion Laboratory (JPL). We fabricate arrays on 350-800 fim 
thick silicon substrates chosen for optimal optical coupling. 
We deposit 0.7-1.2 yum low-stress silicon nitride films (ESN) 
and tune the internal stress to be as low as 150 MPa to ensure 
mechanical integrity. 

The TES normal state resistance Rn and transition tempera¬ 
ture Tc are sensitive to thickness variability and chemical con¬ 
tamination, so we deposit and pattern our TESs before the 
millimeter wave circuitry on a fiat and chemically clean sur¬ 
face. We first e-beam evaporate aluminum (Al) and pattern 
the Tc = 1.2 K TES. We DC-sputter titanium (Ti) and use in¬ 
ductively coupled plasma (ICP) etching to etch the = 0.52 K 
TES for on-sky observing. Both films are immediately chemi- 
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TABLE 1 

Summary of Processing Steps* 


Material/Eunction 

Deposition 

Etchant 

1. Low-stress nitride 

High temp LPCVD 

CHF 3 ICP-RIE 

2. Aluminum TES 

e-beam evaporation 

Liftoff 

3. Si 02 Pro-1 

RF sputtered 

ICP, CHF 3 /O 2 etch 

4. Titanium TES 

sputtered 

ICP, Freon-1202 

5. Si 02 Pro-2 

RF sputtered 

ICP, CHF 3 /O 2 etch 

5. Niobium Ground Plane 

Magnetron sputtered 

Liftoff 

6 . Si02 ILD 

RF sputtered 

ICP CHF 3 /O 2 etch 

7. Gold resistor 

e-beam evaporated 

Liftoff 

8 . Niobium Microstrip 

sputtered 

ICP, Freon-1202 

9. Gold heat capacity 

e-beam evaporated 

Liftoff 

10. Silicon release 

substrate 

STS DRIE/ XeF 2 


* Acronyms defined in text 


cally passivated with RF-sputtered silicon dioxide (Si 02 ) with 
ICP etched via holes to allow DC electrical contacts with 
subsequent layers. The RF bias ensures that the passivation 
patches have rounded edge walls so that subsequent films can 
make contact over steps. The top panel of Figure [7] illustrates 
these first steps. The A1 and Ti films make a series DC con¬ 
nection and are later connected to Nb bias lines. 

The millimeter wave circuits are fabricated from four films: 
a niobium ground plane, a Si 02 interlayer dielectric (ILD), an 
upper microstrip Nb conductor, and thin Au resistive termina¬ 
tion. We DC-sputter the Nb films and pattern the ground plane 
with lift-off, a technique where we first deposit photoresist 
and then dissolve it from under the metal film to be removed. 
In principle, this film could be etched, but we suspect that the 
Ti Tc can be altered by this step despite the intervening protect 
layer. The Nb ground film defines the antenna slots, bandpass 
filter inductors, holes for bolometer release, and safety holes 
under the bond-pads. This film also fills the bias vias down 
to the TESs. Lift-off provides rounded side-walls and thus 
ensures step-coverage of subsequent films. The Nb ground- 
plane covers more than 90% of the 100 mm wafer, serving as 
both the DC and RF ground for the detector array. 

We have explored several ILD materials, although all cur¬ 
rently deployed devices use RF-sputtered silicon dioxide. 
Thickness uniformity of this film is crucial for detector unifor¬ 
mity across the array. A 6 inch Si 02 RF sputtering target and 
substrate rotation help achieve thickness uniformity to bet¬ 
ter than 7% across the array. We use an inductively coupled 
plasma (ICP) reactive ion etch system to etch release holes 
around the bolometers and vias that allow for a DC connec¬ 
tion to the buried TES structures in the subsequent Nb metal¬ 
lization fabrication step. 

We form the resistive terminations on the bolometer islands 
with e-beam deposition and lift-off patterning. Finally, we 
DC-sputter the upper Nb microstrip conductor and etch from 
that film the microstrip feed network, band-defining filter ca¬ 
pacitors, and DC bolometer bias lines and bond-pads using a 
freon etch. We have found that lift-off techniques (used for 
Bicep2) and chlorine-based (BCI 3 ) etches can contaminate 
the Nb films at this step, resulting in unacceptable mill ime- 
ter wave losses in our circuitry as discussed in ^ and ^63 


The middle panel of Figure[7]iliustrates these intermediate mi 
crowave steps that define ground plane, ILD, and microstrip 
traces. 

The final fabrication steps release the bolometer islands, il¬ 
lustrated in The bottom panel of FigurefT] We etch through the 
LSN with CHF 3 in an inductively coupled plasma reactive ion 
etch ICP-RIE system, exposing bare silicon under the release 


viewed above 
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Eig. 7.— Cross-section of films in order of fabrication The aspect ratio 
between radial and normal dimensions are distorted for clarity. We include 
photos of how the device looks face-on for reference. Upper. Deposition and 
etching of films for the TESs and their protect layers. Middle: Deposition and 
etching of antenna and microstrip features. Lower Release of TES bolometer. 

holes. We e-beam evaporate and lift-off thick gold onto the 
islands to add heat capacity and control the required readout 
bandwidth. We use an STS deep reactive ion system to etch 
from the front side and remove silicon from the release holes. 
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This Bosch etching process uses a combination of etching and 
passivation steps to cut vertically through the silicon wafer, 
after which we und ercut the bolometers with a XeF 2 release 
( Turner et al.|2001 ). XeF 2 attacks the remaining exposed sili- 
con isotropically, so the small hole pattern in the photoresist is 
pre-determined to release only the islands and legs and to not 
undercut the remaining antenna structures. Finally, we use the 
STS etcher to cut out the square tile with holes for alignment 
pins. Tablesummarizes these fabrication steps. 

5. MATERIAL PROPERTIES FOR THE MILLIMETER WAVE 
CIRCUITS 


We separately characterize the material properties used in 
the design of the antenna feed and filter, in particular the Nb 
penetration depth and the ILD dielectric constants. These 
material parameters are poorly characterized in the literature 
for millimeter waves, and they can depend on details of the 
processing. We also need to characterize loss properties and 
monitor stability over time. To address these concerns, we 
have fabricated companion tiles of test devices where each de¬ 
tector pair receives power through a single polarization broad¬ 
band integrated antenna. The microstrip feed evenly divides 
power in a micro strip tee-junction between a device under 
test (DUT) and a reference bolometer to divide away opti¬ 
cal effects from the antenna and fore optics. We measure 
the spectral response of the bolometers through Fourier trans¬ 
form spectroscopy(CTS)_mth a Martin-Puplett interferom¬ 
eter ( [Martin & Puplett|[l970| ), combining the detector time- 
streams and encoder readings from the moveable mirror’s 
translation stage to form interferograms. We low-pass fil¬ 
ter, zero-path difference, and Hanning apodize the interfer- 
ogram before Fourier transforming the interferogram into the 
frequency response ^(z^). 


Antenna 


Reference 

Bolometer 


DUT 

Bolometer 



Fig. 8. — Upper. Schematic of a test device for microstrip wave speed, 
where the mismatched microstrip segment at right forms a Fabrey-Perot cav¬ 
ity. Lower. Standing wave pattern in a device’s FTS spectrum. The period is 
proportional to the transmission line wavespeed. 

We use one device for microstrip wave speed measure¬ 


ments. The DUT is a microstrip Fabrey-Perot cavity with 
an impedance intentionally mismatched from the surround¬ 
ing lines. The standing wave pattern in the spectrum of the 
bolometer behind the cavity can be used to determine wave 
speed. (See Figure for a schematic and data sample.) The 
wave speed is a function of both the ILD dielectric constant 
er and the Nb penetration depth A. Alternatively, the band¬ 
defining filter resonance is determined by series inductance 
dominated by the CPW magnetic inductance. Data from the 
band locations and test devices allow us to solve for both 
and A and we summarize these numbers is Table [2] for differ¬ 
ent materials. 


TABLE 2 

Material Properties at millimeter waves 


Parameter 

Value 

Si02 er 

3.9 

SixNy Cr 

7.0 

Nb A 

0.1 p,m 


Another test device measures loss per length in our trans¬ 
mission lines, where the device under test is a stretch of line 
that is several wavelengths long. We determine loss as a func¬ 
tion of frequency by forming the ratio of FTS spectra seen 
through the long line and the reference bolometer. While the 
loss tangent has been shown to be dominated by loss in di¬ 
electrics, and can be quite low in some dielectric systems, our 
measurements determine loss in the microstrip as fabricated. 
Thus far we have tested PECVD Si02, PECVD and 

evaporated Si 02 in this manner. 

This loss test capability became particularly important dur¬ 
ing a period of time when the optical efficiencies of our fabri¬ 
cated devices were unusually poor. This test program allowed 
us to show that the effective loss tangent of our microstrip had 
a frequency dependence (see Figurej^ when using all three 
of the dielectrics above. This common problem present with 
three ILDs from different deposition systems suggested that 
the loss was generated in the Nb films. Ultimately, these mea¬ 
surements helped us to identify and correct a modest tensile 
stress in the Nb as deposited by sputtering. Further, the loss 
appears to be sensitive to the reactive ion etch chemistry used 
in the Nb process. Controlling the film stress at a nominal 
compressive -100 to -300 MPa and changing the Nb etchant 
from BC 13 -F 02 to FI 2 +O 2 resulted in lower loss tangent with 
a dependence (see Figure [^. 

6 . BEAM SYNTHESIS 

The antenna array combines waves from the sub-antennas 
with equal amplitude and phase to synthesize a uniform il¬ 
lumination, which generates a sine pattern in the antenna far 
field. Using a test cryostat, we characterized the antenna far 
field pattern and confirmed that the antenna has a 1 beam 
waist of 4.U (FWHM - lU) and side-lobe level of -12 dB 
relative to peak response, as expected from the design. A sam¬ 
ple far-field pattern is shown in Figure [T^ taken in a cryostat 
without lenses or a stop. For a detailed beams characterization 
in the complete Bicep2 and Keck Array cameras, we refer the 
reader to the Bicep2 and Keck Array Beams Paper (|Bicep2 
[Collaboration IV|2014j ). 

Our ground-based experiments difference detector pairs at 
the time stream level to suppress unpolarized common-mode 
noise from the atmosphere. However, differencing at this 
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Antenna 


Reference 

Bolometer 




Frequency [GHz] 

Fig. 9.— Upper. Schematic of a test device for microstrip loss rates, where 
the DUT is a stretch of transmission line several wavelengths long. Lower. 
Sample spectra showing power loss per unit length (a = -J(logP) /dx) in 
hlms exhibiting low and high loss. Nonlinear increase in loss with frequency 
has been an indicator of bad process parameters in need of adjustment. 
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Fig. 10.— Sample far-held detector pattern measured in a test cryostat 
without an optical stop. Power is normalized to peak on boresight and the 
dashed line indicates where the f/2.2 camera stop would lie. 

early point in the analysis pipeline can allow temperature 
anisotropies to leak into polarization if the detector beam pat¬ 
terns are mismatched. This is an especially acute challenge to 
small aperture experiments because the beam size couples to 
large temperature gradients at degree scales. 

Our team has developed an analysis technique, called de¬ 
projection, that discards contaminated modes. As described 


in the Bicep2 Systematics Paper ( |Bicep2 Collaboration III 
2Q15| ), the pipeline can remove temperature leakage through 
relative-gain mismatch (monopole moment of beams), first 
derivatives of temperature through displaced centroids (dipole 
moment of beams), and second derivatives of temperature 


through differential beam widths and ellipticities (quadrupo- 
lar moments of beams). 

The detector design is most prone to centroid displacement, 
as charac terized in the Beams Paper ([Bicep2 Collaboration 


IV|2014| ). We have identified two mechanisms responsible for 
beam displacement and implemented design and fabrication 
fixes to suppress their contamination below the r = 0.1 level 
( [Q’Brient et al.|[20T^ . These fixed detectors are deployed 
in Keck Array and Spider and Figureshows the centroid 
steer before and after these corrective measures. 

6.1. Parasitic microstrip cross-talk 

The microstrip lines in the feed networks’ horizontal arms 
must fit between the slot sub-radiators, crossing them only at 
the intended feed points. To achieve this, the lines must be in 
close proximity, and for the 150 GHz Bicep2 devices, many 
lines are separated by only ^10 dielectric thicknesses. This 
separation would be adequate to avoid cross-talk for short 
stretches of lines, but some of the line pairs run parallel for 
over four wavelengths. The longest of these pairs are those 
running from the sides to the center to begin the horizontal 
tree arms, coupling power to one of the two lines they inten¬ 
tionally split power between, as seen in Figur e [4 and summa¬ 
rized in the bottom of the left panel of Figure [IT 

In a reverse-time picture, power from port 1 (defined in Fig¬ 
ure [TT]) should be evenly divided between ports 2 and 3. The 
currents in a pair of parallel lines can be expressed as a super¬ 
position of even and odd modes, where the even mode has the 
same current magnitude and directions in each line and the 
odd mode has the same magnitude but opposite directions. 
Coupling between the two parallel lines, however, results in 
the two modes having different field concentrations in the air 
and dielectric volumes, and so induces a difference in wave 
speeds between the modes. Waves passing from port 1 to in¬ 
ternal port (4) thus also couple waves into internal port (3) 
that lag by 90° in phase: 


%1 = ( 5 ^ 4)1 + 5 ^ 4)1 )/2 = 


sm((ke-ko)l/2) 


= -je 


,-Kke+K)tl2 


(9) 

where ke and ko are the even and odd mode wavenumbers and 
£ is the parallel coupling length. This circuit acts as an unin¬ 
tentional reverse-wave coupler, where impedance mismatches 
bet ween even and odd modes are small, but wavespeeds are 
not ( jHammerstad & Jensen||1980| . The waves entering the 
microstrip tee-junction frorn internal port (4) add in quadra¬ 
ture to the intended waves from internal port (3), advancing 
the total phase on the side opposite the vertical tree and re¬ 
tarding it on the other. This phase-step steers the beams off 
boresight away from the vertical summing trees and is where 
the left-right symmetry splitting occurs in the feed. This ef¬ 
fect produces the repeatable horizontal centroid steer seen in 
the cameras’ near field. We have reproduced this in simula¬ 
tions using HFSS that account for the upper conductor finite 
thickness. 

As seen in the right panel of Figure[TT] we have re-designed 
the antenna feeds for the Keck Array and Spider detectors 
to greatly increase spacing between the lines and therefore 
reduce the parasitic coupling effect. We also include ad¬ 
justable phase-lag lines before each slot to remove residual 
phase-error and to synthesize matched beams, although the 
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Eig. 11.— Details of microstrip layout in the center of the antenna feed. Left: Center of the BiCEP2-era design. The Cartoon under the picture shows the 
effective coupler circuit and the ports corresponding to the scattering parameters in Equation]^ The arrows are phasors that illustrate how the intended waves 
(horizontal) combine with the parasitic coupled waves (vertical). Right: Modern design with r^uced coupling due to greater spaced lines. The extra lengths of 
line that form compensating phase lags are visible (e.g. bottom of picture), although these are not necessary with the suppressed microstrip coupling. 


current antenna-feeds’ increased spacing renders the correc¬ 
tive phases unnecessary. These phase-lags are visible in the 
bottom of the right picture in Figure pT] 

6.2. Niobium contamination 

Magnetic fields can penetrate into a superconductor by a 
characteristic depth Impurities in niobium films scatter 
Cooper pairs with a mean free path increasing the pene¬ 
tration depth beyond the London depth Al ~ 50 nm of pure 
niobium to \eff = Al where ^ 40 nm is the Cooper 

pairs coherence length. As displayed in Table our films are 
typically measured to have Xeff ^100 nm, suggesting mean- 
free path of ^ = 10 nm, well within the “dirty” limit where film 
cleanliness can impact circuit performance. 

Nonuniform contamination can produce nonuniform ki¬ 
netic inductance, which can spatially perturb the wavespeeds 
in the microstrip summing tree. Variations in wavespeed can 
steer beams off boresight. Additionally, the summing tree 
does not treat the two polarizations identically, and as a re¬ 
sult, they can be steered differentially. We have observed that 
the tiles with the largest scatter in beam centroid position cor¬ 
respond to those with the largest vertical dipole components. 
We also expect the tree to induce larger steering in the verti¬ 
cal than horizontal because slots along rows combine immedi¬ 
ately in the horizontal tree, resulting in less integrated phase 
error than those along columns that combine in the vertical 
tree after horizontal summing. 

We have built models for how our antenna-feeds perform 
with film gradients. We subdivide the circuit into short sec¬ 
tions of transmission lines and tee-junctions and use each sec¬ 
tion’s location in the detector to assign unique film properties. 
We construct simple scattering matrices for each section and 
casc ade them into one large 289 x 289 matrix per polariza- 
tion( |Bodharamik et al.|1971| ). From this we can compute slot 
illuniiination patterns and thus far-field patterns. We find that 
20% variations in Xeft across the array can produce differen¬ 
tial pointing that is 10% of the beam FWHM in the vertical 
directions, and half that on the horizontal, similarly matching 
our observed scattering in pointing. 

We defined our microstrip lines in early tiles with the lift-off 
technique that we use for the Nb ground plane, as described 


in ^ Several devices have shown discoloration in this step, 
leading us to speculate that the Nb leaches organic materials 
from the resist during lift-off. These observations and model¬ 
ing inspired a switch to an etch-based means of defining the 
Nb microstrip lines (described in Q. This simple fix reduced 
the scatter in centroid location to ^ 1% of Gaussian width a, 
and the right panel of Figure shows the centroid alignment 
between polarization pairs foFboth colors of deployed focal 
planes. 

6.3. Direct stimulation of bolometers 

Optical power is meant to reach the bolometers only 
through the antenna and microstrip feed network. However, 
photons can directly excite responses in the detectors which 
is of particular concern in a design where there are no horn 
blocks to shield the sensors themselves. In Bicep2 pre¬ 
deployment testing of early generations of detectors, response 
to out-of-band power was detected at levels 3%-4% of the 
total response, with near field and far field angular response 
patterns consistent with direct stimulation of the bolometer is¬ 
lands. Simulations in HFSS and CST suggest that the ground 
plane on the bolometer islands can be inductively held at a dif¬ 
ferent voltage than the surrounding ground plane through the 
microstrip ground that we deposited only on one leg in early 
generation prototypes. This voltage can drive millimeter- 
wave currents through the gold-termination, resulting in a di¬ 
rect stimulation of the detectors not through intended antenna 
and microstrip feed. Steps taken prior to Bicep2 deployment 
to minimize this coupling included the addition of metal mesh 
low-pass edge filter above the focal plane and several design 
changes to the bolometer islands themselves. The island leg 
design was modified to narrow the width of the opening in the 
ground plane surrounding the island, and ground plane conti¬ 
nuity was extended onto the island by metalization of the four 
outer support legs. 

With these modifications in place, we have taken Bicep2 
and Keck Array 150 GHz camera optical efficiency measure¬ 
ments with and without high-pass “thick-grill” filters that ob¬ 
struct power below 200 GHz; after accounting for the filter’s 
filling factor, we found that the optical response through the 
filter was ^ 0.5% of that without filtering. This small leak- 
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Horizontal Offset [inches] Horizontal Offset [arcmin] 


Fig. 12.— Left: Beam centroid displacements in the near field before a nd after fixes were mad e to the etch recipe. The performance sh own in left panel would 
allow a BiCEPl-style pipeline without deprojection to constrain r to 0.1 JTakahashi et al.|2010|. The deprojection pipeline described in JBicep 2 Collaboration! 
|ni|2015) allows r to be constrained to yet lower levels. Right: Beam centroid displacements m the far field for two different colors after the recipe was fixed. 


age suggests limited “blue-leaking” that would allow above¬ 
band power to excite response in the bolometer. We have 
also measured the response of our detectors to a chopped ther¬ 
mal source on boresight through a polarizing grid at different 
angles; we have found that the crossed response is similarly 
^ 0.5% of the co-polarized response, consistent with known 
multiplexer cross-talk levels. If the detectors are acting as 
a direct absorber, then their small area compared to the an¬ 
tenna’s should provide a broad angular response. As a result, 
we expect that the //2.2 stop in the cameras for these exper¬ 
iments helps limit the direct stimulation and we have found 
that tests of devices in cameras with faster optics can have 
higher direct stimulation levels. 

7. ARRAY PROPERTIES AND UNIFORMITY 

Our high attained sensitivities are a result of our high detec¬ 
tor yield, a figure that depends on attaining uniform detector 
properties across the array. This uniformity also helps to mit¬ 
igate some potential sources of systematic error. This section 
describes array uniformity of a variety of properties. 

7.1. Bolometer thermal conductance G 

Thermal fluctuation noise across the bolometer legs is the 
largest internal source of noise in Bicep 2, and is given by: 

NEPg" = AkT^GcF{Tc, T„) (10) 

where the “Mather factor” FiJc^To) accounts for thermal gra- 
dients across t he bolometer legs and varies between 0.5 and 1 
( Mather] 1982 1. The saturation power of the detector is given 
in Equation [8 Note that the leg thermal conductance G is a 
function of temperature. 

The legs’ total parallel thermal conductance G needs to be 
chosen to keep thermal carrier noise subdominant to photon 
shot and Bose noise, but avoid saturation under on-sky load¬ 
ing. We try to have the total saturation power (Joule heating 
from the bias circuit plus optical loading on the detectors) ex¬ 
ceed the optical loading (^4-6 pW at 150 GHz in Bicep2) 


by a factor of two. This aggressive safety factor is only possi¬ 
ble if we have tight control of our processing parameters that 
determine G and Tc (discussed in the next subsection). Our 
achieved repeatability between detectors in a tile and from 
tile-to-tile are shown in Figure Our time-division multi¬ 
plexing readout requires that detector s in a common readout 
column be commonly voltage biased (|Battistelli et al.| (|2008 ) 
and [Stiehl et ak] ( |2Q11| )), and our control of G within a tile 
allows this. 



[pW/K] 

Fig. 13.— Thermal conductance at Tc. This figure demonstrates repeata¬ 
bility within and between tiles. 


The XeF 2 release only has a 10:1 selection ratio between 
silicon nitride and silicon, and variation in bolometer leg 
thicknesses from nitride etching would seriously compromise 
these numbers. For this reason, we oxidize the wafers prior 
to nitride deposition. The thin ^100 nm Si02 film has a 
much higher 100:1 selection ratio between oxide and sili¬ 
con and th us protects the legs from XeF 2 attack ( [Williams 
et al.|2003] ). This oxide film ensures more uniform leg cross- 
sections across the detector arrays, and thus more uniform 
bolometer Gs. 
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11. Transition temperature Tc 

The balance between the between the internal ^He/^He/'^He 
fridge cooling power and thermal loading lets the focal plane 
cool to ^ 250 mK. We use heaters to operate the focal plane 
at an elevated 280 mK so we have margin to control this tem¬ 
perature through active feedback loops while still maintaining 
photon-noise-limited sensitivity. This choice provides the de¬ 
tectors with a bath temperature Ti, = 280 mK. For bolometer 
legs with n=2.5, NEPg is minimized at an island tempera¬ 
ture of 420 mK, although this is a broad-optimum with a gen¬ 
tle slope at higher temperatures. For on-sky operation, elec¬ 
trothermal feedback locks the bolometer to the titanium TES’s 
Tc ^ 520 mK. While this temperature is higher than optimum, 
it provides some extra margin against our detectors latching 
into a pure superconductive state while only increasing the 
NEP from thermal conduction by 4%. By selec ting a pure su 


perconductor for this TES instead of a bilayer ([Fabrega et^ 
( [2009| ), |Gildemeister et al.| ( 19991 ) and JMyers et al.| ( |2005[ )), 
we obtain repeatable transitions devoid of the multiple tran- 
sitions that can result from a bilayer’s complicated chemistry. 
Deposition on bare wafers and passivation of oxide protection 
layers also help maintain uniform values between and within 
tiles. As seen in Figureour titanium Tc is uniform across 
tiles and even between tHes, maintaining uniform saturation 
powers (Equation within a multiplexer column, necessary 
to find a common bias point. 



7.3. Normal resistance 

TES bolometers will only experience electrothermal feed¬ 
back if voltage biased, and we bias our nominally ^ 60 mfl Ti 
TESs with a 3 mQ shunt resistor. The Al TES is even higher 
in resistance, so both are far larger than the shunt. However, 
this condition is only met if the Ti TES resistance is repeatably 
this large. We also need them repeatable in value so the detec¬ 
tors will have a common latch resistance and thus a common 
bias point where they will experience strong feedback. Fig¬ 
ure El shows that these values are indeed uniform between 
and within tiles. Just as was true for Tc uniformity, our recipe 
of a pure TES deposited as first layers on the tile, protected 
with an oxide film, helps maintain repeatable performance. 


7.4. Time constant r and loop gain ^ 

The effective thermal time constant of an ideal voltage- 
biased TES bolometer is given by 


G/C 

\+^{vy 


( 11 ) 



Fig. 15.— Titanium Normal Resistance for four Bicep2 tiles 
where the intrinsic bolometer time-constant C/G is decreased 
by the effective loop gain of electrothermal feedback, The 
responsivity of our detectors is the change in measured current 
for a change in incident optical power, given by 


-1 ^ 1 


( 12 ) 


Both of these expressions experience corrections for non¬ 
idealities such as finite shunt resistance. Equations 0 and 


12 indicate that a high loop gain ^ will increase the speed of 
the sensor and simplify its responsivity (s ^ -1 /Vt). Operat¬ 
ing the bolometers in this limit also maintains a fast detector 
response (low r), giving rise to an approximately fiat detec¬ 
tor transfer function for the frequencies of int erest, which for 


our ground-based experiments are / < 2 Hz (Irwin & Hilton 
|2005l ) 


In order to maintain stability against electrothermal oscilla¬ 
tions, the detectors’ thermal bandwidth must not exceed that 
of the electrical bias circuit. The Bicep2 bias circuit includes 
a 1.35 pH inductor in series with the ^ 3 mfl bias resistor to 
avoid aliasing above a roll-off of R/L 5-6 kHz. To limit 
the TES bandwidth, we deposit thick ^ 2 pm gold onto the 
Bicep2 detectors that adds an additional C ^ 0.5 pJ/K heat 
capacity to the bolometer island. The resulting time constants 
at typical science biases are shown in Eigure jl6| 

Spider detectors are intrinsically slower due to their lower 
leg conductance. It is thus not necessary to add as much gold 
to the island to ensure stability; typically only ^0.5 pm is 
deposited. Time constants on transition are found to be sim¬ 
ilar to Bicep2’s. Eor Spider detectors we have measured 
time constants in response to optical square-wave excitations 
at a range of bias voltages (and thus TES resistances), allow¬ 
ing us to infer that the normal time constants (without loop 
gain) are G/C ^ 30 ms and that loop gains are ^ ^20-30 
Sit R = 0.6Rn (Eigure[T^. Bicep2 TESs should have similar 
loop gains. We note m passing that fast bolometer time con¬ 
stants are useful for limiting the effect of particle radiation 
(e.g., cosmic rays) on detector time streams in balloon- and 
space-borne instruments. 


7.5. Spectral Response 

Our cameras use a series of low-pass filters, both absorp¬ 
tive plastic and reflective metal-mesh, to limit thermal loading 
on the focal plane and above-band response of our detectors. 
However, we rel y upo n the integrated microstrip filters de¬ 
scribed in Section [3^ to avoid the atmospheric lines immedi¬ 
ately adjacent to our observing bands. Eigurep^ shows mea¬ 
sured response ^(z^) averaged across a focal ^ne for three 
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Fig. 16.— Upper: Time constants of Bicep2. These were only carefully 
measured in transition. Lower: Loop gain of Spider. We measured time 
constant at a variety of biases for Spider characterization, allowing a proper 
measurement of loop gain. Both sets of detectors had similar r ~ 0.9 ms 
biased time-constants. 



Fig. 17.— Measured detector spectra for devices designed for 95 GHz 
(red), 150 GHz (green) and 220 GHz (blue). The data for the 150 GHz are 
from Bicep2, while the others are from a test cryostat. Winter atmospheric 
transmission at the South Pole is overlaid in black. 

different spectral channels (in three different focal planes), 
demonstrating that this technology does indeed avoid atmo¬ 
spheric features. 

We can summarize the spectral response ^(z^) by comparing 
it to a top-hat response with the same area and then defining 
the band-center as: 



h'S(a)diy 


and the bandwidth as 


^ f sHv)dv ■ 


(13) 


(14) 



Fig. 18.— Spectral response of detectors in the Bicep2 camera. Upper: 
histogram of the band centers. Middle: histogram of band widths. Lower: 
histogram of absolute calibration (abscal) mismatch, where difference and 
mean are computed between polarization pairs. Calibration is through cross¬ 
correlation against Plank 143 GHz. 


The spectral responses of the detectors within a polariza¬ 
tion pair must be closely matched to avoid systematic con¬ 
tamination through spectral gain mismatch. In principle, 
using separate filters for each polarization might induce a 
mismatch through gradients in material properties. Our de¬ 
sign resists this by placing the filter pairs in close physical 
proximity; in the 150 GHz detectors, the filters are mm 
apart. Furthermore, the filters’ magnetic inductance domi¬ 
nates the kinetic inductance; this renders our filters’ features 
particularly robust against variation in Nb contaminati on that 
might adve rsely impact electrically long resonators ( [Myers | 
et al.||2005] ). As shown in Figure our band centers and 


widths are found to be highly repeatable. We calibrate our 
responsivity through correlation against the Planck 143 GHz 
maps. Spectral mismatches between filter pairs would mani- 
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fast themselves as differences in each detector’s absolute cali¬ 
bration (abscal), but the bottom panel of Figure shows that 
these are less than 3% of the average responsivity, with a me¬ 
dian of 0.5%. 


7.6. Optical Efficiency 

We characterize our end-to-end optical efficiency by bias¬ 
ing onto the aluminum TES, which has a higher Tc than Ti, 
and thus higher saturation power. By comparing the measured 
response to aperture-filling 77 K and 300 K sources to the ex¬ 
pected incident power Fine = kAT f f(a)diy, we can use the 
measured spectra to infer efficiency. We routinely obtain total 
camera optical efficiencies in excess of 30%. For example, 
the Bicep2 efficiencies are shown in Figure [T^ These fig¬ 
ures are for the total camera efficiency, which includes losses 
other than the detectors. Measurements of response of early 
engineering-grade detectors to an internal cold load, which 
exclude spillover loss and loss in our filter stack, su ggest that 
the raw detector efficiency should be nearly 70% dOrlando 
let al.|2009D . 



Fig. 19. — Optical Efficiencies of Bicep2 detectors. Blue curves (left-axis) 
are end-to-end receiver efficiency through all optics; green curves (right-axis) 
are raw detector efficiencies for a single test-tile from an engineering-grad 
test focal plane, in response to an internal cold-load. 


8 . SENSITIVITY 

The array design described above was first deployed for 
astronomical measurements in 2009 as part of Bicep2, a 
ground-based CMB polarimeter. We have since fabricated 
dozens of deployment-grade arrays at 95 and 150 GHz for 
use in the terrestrial Keck Array and balloon-borne Spider in¬ 
struments. These programs have yielded extensive data on the 
real-world performance of this technology, as well as demon¬ 
strating its adaptability to different optical loads. 


8.1. Noise 

The th eory of noise in TES bo lometers is by now well de¬ 
veloped prwin & Hilton|[?005| ). Major contributions typi¬ 
cally include photon noise, phonon noise (thermal fluctuation 
noise), Johnson noise from the TES and shunt resistances, and 
amplifier noise. In many cases TESs also exhibit varying de¬ 
grees of “excess” noise beyond that predicted by simplified 
models. In multiplexed systems, noise performance also de¬ 
pends upon the relationship between the noise levels and de¬ 
tector and readout bandwi dths, since poor choices can lead to 
substantial n oise aliasing ( IBattistelli et al.| ( 2008| ), and |Stiehi 
|et al.| ( |2011| )). We find that the detectors’ measured noise 
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Fig. 20.— Measured and modeled noise for a single detector in the Keck 
Army. The red line indicates the Nyquist frequency for the multiplexing rate. 


is well reproduced by a simple model incorporating mod¬ 
est excess TES noise, and we have successfully operated the 
bolometers in configurations with little aliasing penalty. 

In order to characterize the noise performance, we have 
measured noise spectra for the Bicep2 and Keck Array de¬ 
tectors using un-multiplexed data digitized at 400 kHz. Al¬ 
though the frequency range corresponding to degree-scale 
anisotropies is 0.05-1 Hz for the Bicep2 and Keck Array scan 
strategy, these measurements allow us to observe device per¬ 
formance near the typical multiplexing frequency of 25 kHz 
and thus model noise aliasing due to the readout. These mea¬ 
surements, adjusted for expected aliasing, are in good agree¬ 
ment with noise measured in the science-mode multiplexer 
configuration. 

Figure ^compares the measured noise for a representative 
detector from the Keck Array to its various modeled compo¬ 
nents. Photon noise dominates at the low frequencies of inter¬ 
est for sky observations. The photon noise equivalent power 
(NEP) can be expressed as a sum of Bose and shot noise con¬ 
tributions: 

20 ^ 

NEP^hoton = + (15) 

where u is the band center, Az/ is the bandwidth, and Qioad 
is the optical loading. The next-largest contribution to the 
noise at low frequencies is the thermal fluctuation (phonon) 
noise across the SiN isolation legs, given by Equation [T^ All 
other modeled contributions, including the TES Johnson and 
excess noise, the shunt resistor noise, and the cold and warm 
amplifiers, are negligible at low frequencies. 

At frequencies of « 1 kHz the TES excess noise starts 
to contrib ute significantly. The T ES excess noise, as de¬ 
scribed in ( |Gildemeister et al.|2001| ), tends to be proportional 
to the TES transition steepness ^ {dl/dR) |r, which 

for these detectors is higher at lower resistances. The excess 
noise varies between fabrication batches, although the detec¬ 
tor in Eigure j20| has a relatively low amount. The excess noise 
should have httle contribution to an experiment’s overall noise 
level as long as it is not aliased. Bicep2 and the Keck Array 
avoid such aliasing by multiplexing at 25 kHz and biasing 
at relatively high resistances on detectors that have large ex¬ 
cess noise components. Spider cannot multiplex as quickly 
due to long cable lengths, so we use a lower TES resistance 
{Rn ^ 30 mft) and higher Nyquist inductance (2.0 pH) to limit 
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TABLE 3 

Summary of detector noise are low erequencies eor Bicep2 and Keck Array. 


Receiver 

Bicep2 

Keck Array RxO 

Rxl 

Rx2 

Rx3 

Rx4 

Absorbed Photon Noise [aW/ \Tlz] 

41 

33 

32 

34 

27 

27 

Phonon Noise [aW/v^] 

27 

24 

20 

24 

14 

16 

TES Johnson Noise [aW/v^Hz] 

0.5 

0.8 

0.5 

0.8 

0.4 

0.6 

Amplifier Noise [aW/v^Hz] 

2.0 

2.0 

2.0 

2.4 

2.0 

2.5 

Total Noise + Aliased at 25 kHz [aW /\/Hz] 

56 

46 

39 

48 

36 

40 

per-detector NET from maps [/xKcmb V^] 

305 

316 

351 

317 

370 

429 


aliasing of excess noise. 

The noise contributions for Bicep2 and the Keck Array at 
low frequencies and under winter atmospheric loading condi¬ 
tions are broken down in Tablefor the observing year 2012. 
The photon noise is of absorb^ power only to aid the com¬ 
parison between models and measurements without adding an 
extra factor of optical efficiency. The last line of the table 
shows total NET computed from jackknife maps, which nec¬ 
essarily references noise of the incident photons and accounts 
for the optical efficiency of the cameras. 

8.2. Performance of ground-based designs 

We have measured the devices’ performance on the sky dur¬ 
ing lengthy Bicep2 and Keck Array observation campaigns. 
The sensitivity of the experiment can be measured using the 
timestream noise between 0.1-1 Hz, after calibrating to the 
CMB. Bicep2 had a noise equivalent temperature (NET) of 
16 /iKcMB\/s- The larger detector complement of the Keck 
Array achieved 11.5 and 9 /iKcMB in 2012 and 2013, re¬ 
spectively, also at 150 GHz. Eigure^ shows the distribution 
of per-detector sensitivity for the KeHt Array in 2013. 


450 



400 600 

jlK y/s 


1000 


Fig. 21.— NET per detector histogram for the Keck Array in 2013 


Another measure of instrument performance is the total 
map depth achieved, defined as the noise within a map pixel 
of a specific size. This is a function of instrumental sensitivity, 
scan strategy, and observing time. Over three seasons of ob¬ 
serving, Bicep2 achi eved an rms noise of NEQ = 87 nK - deg 
(5.2 //K • arcmin). (|Bicep2 Collaboration iT] ( |2014| ), |Bi-| 
CEp2 Collaboration I| ( [2014] ))! Over two years, the five 


Keck Array 150 GHz cameras achieved NEQ = 74 nK • deg 
(4.4 yuK • arcmin). Combining Bicep2 and Keck Array and 
averaging across the entire effective area 400 square degree 
field results in rms noise of 2 nK. 

9. CONCLUSIONS AND FUTURE WORK 

We have demonstrated our novel detectors’ performance 
through 18 camera-years of observations in the BlCEP2/Kpck 
Array program, including a detect ion reported in 2014 
of degree-scale B-mode anisotropy (|Bicep 2 Collaboration 
|il W- This paper has described the design princi¬ 
ples, challenges, fabrication techniques and our characteriza¬ 
tion/screening program that made these deployments possi¬ 
ble. This program has continued through the recent, success¬ 
ful 2015 Spider flight and higher frequency upgrades to Keck 
Array. 

We have also deployed Bicep3 in the current (2014-2015) 
Antarctic summer season, which has 1152 95 GHz detectors 
in a single camera. Once all the detectors are installed in the 
2015-2016 season, this instrument will support 2560 95 GHz 
detectors. We package these detectors in individual modules 
that efficiently fill the focal plane and we illuminate the de¬ 
tectors’ antennas with non-uniform Gaussian tapered slot il¬ 
luminations to reduce spillover onto the camera’s stop. These 
recent modifications will be the subject of a future paper. 

Lastly, our detector technology allows for multi-color focal 
planes where different color channels are co-located on the 
focal plane. Most competing dual-band detector technologies 
use a single common aperture for all color channels and thus 
there is a reduction in per-dete ctor efficiency due to aperture 
spillover ( [Q’Brient et al.|201^ . Planar antenna arrays allow 
each color to have a custom aperture, thus more efficiently us¬ 
ing of both focal plane real estate and limited readout capac¬ 
ity. These are under active development for Bicep3 further 
upgrades to the BlCEF/Keck Array program, and for future 
Spider flights. 
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